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Abstract 



The dark matter accretion theory (around a central body) of the author 
on the basis of his 5-dimensionaI Projective Unified Field Theory (PUFT) is 
applied to the orbital motion of stars around the center of the Galaxy. The 
^ ' departure of the motion from Newtonian mechanics leads to approximately 

, flat rotation curves being in rough accordance with the empirical facts. The 

' spirality of the motion is investigated. 

^ . 

2 ■ 1 Retrospect 



In a preceding paper (Schmutzer 2001) we developed a theory of dark matter accre- 
tion around a central spherically symmetric body on the basis of PUFT (Schmutzer 
1995, 2000a, 2000b, 2000c) which we applied to the motion of a test body around 
a central body with the result of an additional radial acceleration effect towards 
the center, compared with the motion on the basis of Newtonian mechanics. Par- 
ticularly this effect was used to explain the empirically measured deviation of the 
. motion of the satellites Mariner 10/11 (Pioneer effect) from the Newtonian motion. 

^ ' In this paper we apply our dark matter accretion theory to the orbital motion 

of the stars around the center of the Galaxy with the aim to find an explanation of 
the approximetely flat rotation curves of the orbiting stars. 

In this context we mention interesting papers of Dehnen et al. (Dehnen, Rose 
and Amer 1995) who elaborated a theory of the flat rotation curves on a different 
theoretical basis applying the Bose-Einstein statistics (preferred by physical argu- 
ments) and the Fermi-Dirac statistics (dropped by physical arguments). By these 
authors quantum-statistical calculations were performed in great detail. 

For a rough understanding of our approach we restricted our investigations to the 
Boltzmann-Maxwell statistics describing a gas of one sort of dark matter particles 
(dm-particles) without degeneracy. We think this step is legitimate to simplify the 
sitation, with the advantage to be able to perform the calculations fully analytically. 

According to this dark matter accretion theory, applied to a spherical homoge- 
neous central body following relationship between the gravitational potential x(r) 
and the dm-particle number density n(r) reads: 

nmx , , 

(n dm-particle number density at "infinity" (great distance), m mass of a dm- 
particle, T kinetic temperature of the gas, k Boltzmann constant). 
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The potential obeys the differential equation {jn Newtonian gravitational con- 
stant) 

Ax + = 47r7Ar/i , (2) 

where with respect to the mass density the quantity /t means /zq (mass density in 
the interior) and (mass density in the exterior). Further k'^ is defined by {x 
dark matter parameter) 



2 47r7ArTO^n 

(3) 

One should realize the positive sign in front of the second term on the left hand 
side of equation (^. Let us in this context mention that we arrived at this approach 
to our theory of dark matter accretion by our experience in the Debye-Milner theory 
of strong electrolytes, where despite attractive forces between negative and positive 
ions at this place a minus sign appears which leads to a Yukawa-type potential, 
whereas in our theory periodic functions occur. 

For presenting the gravitational potentials etc. in the interior and in the exterior 
of a homogeneous sphere considered following two basic functions are useful: 

a) U{z) — z sin z -f- cos z , 

(4) 

b) V{z) = sinz — zcosz 

which determine the potentials as it will be shown in the following. 

For the gravitational potential and the corresponding derivatives we received for 
the interior (index i) and the exterior (index e) of the sphere the results (tq radius 
of the sphere) 

a) Xi ^ Bo{l - —U{}<ro)sm{Kr)} , 
xr 

(5) 

ar >fr^ 

a) Xe = -—V{M:ro)cos{M:r), 

(6) 

b) % = ^v{xro)U{xr). 

dr Hr'^ 

The constant of integration Bq occurring in these formulas is given by 

-DO - — — - — , [/) 
m n k-^Vq 

where 

Mc = ^ (8) 

{Mc mass of the central body) holds. 

Let us for application in the following mention the equation of motion of a test 
body 

dv dlncTp 

- + gradx + .^ = 0, (9) 

where v is the velocity of the test body and (Tc is the cosmological scalaric world 
function. 
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2 Motion of a test body in polar coordinates 

Now we use the polar coordinatets {R, $ }. Then after integration the equation 
(||) reads (dot means derivative with respect to the time t) 

a) R-R^'^ + ^ + RT.^O, 

di? 

b) i?^$(Tc X v\<Jc = Fo with (10) 

„ dinar 
^) ^ = 

As it is well known, the azimuthal velocity is defined by 

z;$ = i?$=-^. (11) 

R(7c 

With the aid of this expression formula (p^a) leads to 
Hence the equation 



v^ = \IRR + RRY, + R^ (13) 
V dR 

follows. Now we introduce the abbreviation 

yl = i? + 7?E. (14) 

In the usual theory of the rotation curves on the basis of Newtonian mechanics 
approximately the motion of the body on a circle (R = const) is considered. We 
think this approximation is not good enough. Therefore, induced by the empirical 
situation for small values of the radial coordinate R we try the basical ansatz 

A = aR^ 10-^"aR, (15) 

where a or a (the factor 10^"^° is an adaption factor with respect to the large 
numerical astronomical values) is a quantity to be specified by further concrete 
physical assumptions . Substituting this expression into (^3|) gives the formula for 
the so-called rotation curve V — w$(i?): 



Finishing this section we would like to emphasize that in our theory the situation is 
basically different to the usual theory of the rotation curves based on the assumption 
R — const, since here because of the periodicity occurring in the potential the 
dx . . . . 

expression - — is indefinite, i.e. that for A = the radicand exhibits negativ-valued 
dR 

regions. 
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3 Rotation curve for the interior and exterior of 
the central body in general (azimuthal motion) 

Interior (index i): 

With the help of (^) and equation (|l^) takes the shape (approximation by 
power series expansion with respect to (>fi?) <C 1) 



=R^a,+ ^^'^^^3^3 U{^ro)V{>cR) . (17) 
Hence for the angular velocity (differential rotation) follows 



Exterior (index e): 

In this case with the help of (^j) and (0) equation (|l6|) takes the shape 



V, = Rja, + V{xT^)U{xR) . (19) 

Analogously follows here the angular velocity 



4 Physical processes in the interior of the central 
body (radial motion) 



Our next task is to get some information about the quantity A in the relation (|16|) by 
investigating physical main processes in the interior of the central body. We propose 
to approach this problem tentatively by considering two opposite processes: 

1. Expansion process of the central body by scalaric cosmological bremshcat 
production according to our hypothesis (Schmutzer 2000c). 

2. Scalaric- adiabatic cosmological contraction process according to our approx- 
imate treatment of the orbital motion of a test body considered (Schmutzer 
2000a). 

4.1 Expansion process (radial motion) 

In the paper quoted above we started from the radial expansion formula of a sphere 
(r radius of the sphere) 

drg = iacrdT (21) 

(T temperature of the body, Oic cubic heat expansion coefficient, index Q refers to 
heat). Further the relationship 

dg = WdT (22) 
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holds {Q produced heat), where for the heat capacity 



W = cqM 



(23) 



(M mass of the central body considered here, eg specific heat of the body) is valid. 

Let us now tentatively identify the heat Q with the scalaric cosmological bremsheat 
Qb produced by the bremsheat effect predicted by us: 



dg = fsdQs 



(24) 



(/s free scalaric consumption factor). Two years ago we derived for the differential 
of the bremsheat production the expression 



dQs = dQiJ™"^')+dg(,™*) 



a' 



din (Tc 



(25) 



(M mass of a body with translational and rotational motion, rc radius vector to 
the center of mass of the body considered, lat moment of inertia tensor, uja angular 
velocity, indices a,b — 1,2,3). For a massive sphere rotating around a fixed axis 
results 



dQB = 



Mri 



2M 



din (Jc 



Using the abbreviation 

J = Mrl; + —{rou;f 
we can simplify (p6|): 

dg B = Jd In CTc . 
Then by means of (^) from ( p2[ ) results 

dT=:^dlna,. 

Now it is convenient to introduce the scalaric-thermal velocity coefficient 



3W 



(26) 



(27) 



(28) 



(29) 



(30) 



Then with the help of (|2^) we receive from (|T|) the bremsheat induced radial 
velocity 



dt 



aqr . 



(31) 



Here one should remember that the coefEcent aq as well as other cosmological 
quantities exhibit a slight scalaric-adiabatic cosmological time dependence. 
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4.2 Scalaric-adiabatic contraction process (radial motion) 

As quoted above, some years ago we treated the circle-Uke orbital motion of a test 
body, arriving on the basis of our scalaric-adiabatic approximation at a contraction 
process with following result for the orbital radius of the moving body: 

= A. ^! 2 n • (32) 
Hence we get by differentiation the expression 

We would like to remind the reader that according to PUFT the true constant mass 
of a body is the (invariant) mass Ai (calligraphic) which is related to the (inertal) 
mass M and to the (gravitative-scalaric) mass M*^'^^ as follows (here a is the general 
scalaric field function): 

a) M = Ma, 

b) M(^) = A^cr(^) with 



(34) 



c) cr^'^-' = — ^ (scalaric force correction factor). 



As pointed out in our publications, since about 1000 years after the big start the 
inequality cr^ ^ 1 is well justified, i.e. the equations (|3^) and (33) read much 
simpler: 

= ff 3 (35) 

and 

rs- ir-,-^- (36) 

By differentiation we find the expression 

rs = ^ ( - -a!) . (37) 



^nMcCTc 

Substituting ( ^5|) till ( |37| ) into the relation ([T^), being applied to this case, yields 
for the rearranged formula (jlj) the result 

As 3 /. 3 .2\ 

as = — = o-c o-^ • (38) 

rs o-c \ CTc / 

For eliminating the second order derivative of the scalaric field function ac we 
have to go back to the corresponding field equation contained in our publications 
(Schmutzer 2000b), namely {L rescaled radius K of the cosmos, rj rescaled time t, 
1? energy density of the substrate in the cosmos, Aq = lO^^cm rescaling factor, xq 
Einsteinian gravitational constant) 

N 3 dL dcTc . , , 

a) „ — 6XV — - — — with 

d?7^ L aq &q 

(39) 

b) i<= = 3.4617 -10 V'cms'. 
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In this context one should remember the rescaling relation 
After this recalculation the formula reads 



(40) 



as 



1 / duc 
dry 



1 dL d(Tc 
L drj dry 



(41) 



Finishing this subsection we eliminate Fq in (pq) by means of (p5|) and obtain 



drc 



-3rd In ctc 



(42) 



4.3 Superposition of both radial effects 

In both preceding subsections we treated in detail the expansion effect ( ^ ) and the 
contraction effect (|4|). By adding both effects we arrive at the total effect for the 
interior 

dr, = dvQ + drs . (43) 
Introducing the abbreviation 

5 = ag - 3S (44) 
leads us to the relationship between the radius r and the total radial velocity 

r,; = ari . (45) 
For further calculations it is convenient to use the scalaric bremsheat parameter 
J fsac 



2,W ' 

Then ^ and (||) take the form 
a) aq — PqT, and 



(46) 



(47) 



b) 5 = (/3q-3)S. 

For the following we specialize our considerations, having led to relation (p6|), to 
the central body (M — > Mc) and further particularly to the case of a central body 
with a resting center of mass {rc'= 0). Then formula ( p7|) reduces to 



J 



(48) 



Neglecting the scalaric-adiabatic time dependence of the material constitution quan- 
tity Pq , then by means of (pTpp) we arrive by differentiation of ( ^ ) at the relation 



= (/3q-3) I] + (/3q-3)e2 



(49) 



Inserting the expressions (^^ and (|9|) into the formula (|i^), applied to the 
interior, leads to (i? ^ r^) 



A, 



air I 



with 



(/3q-3) S + (/3q-2)I]^ 



(50) 
(51) 
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5 Physical processes in the exterior of the central 
body (radial motion) 

The difference between tlie exterior and the interior consists in the fact that in the 
exterior no bremsheat effect takes place. This means that a moving body firstly, 
when leaving by jet the central body, starts with the constant initial velocity 



r-jet = rQ{r = ro) = ag , 

and secondly underlies the scalaric-adiabatic cosmological contraction 
position of both effects gives 



In this case according to (14) follows 

Ae = fe + = OL^T , 

where 



with 



A 



2^0 



(52) 
Super- 

(53) 
(54) 

(55) 

(56) 



holds. 



6 Orbital motion of a test body in the exterior of 
the central body 

6.1 Rearrangement of the velocity equations in usual astro- 
nomical units 

In presenting the plots of the rotation curve and the spiral motion, in astrophysics 
the azimuthal velocity of the moving body (in astrophysics the notion "rotation 

velocity" is used) is indicated in , whereas its distance from the central body is 

given in kpc. In adapting to this practice we have to rearrange the equations ( |l9[ ) 
and (p3|), using the relationship Ikpc = 3.086 • 10^^ km 
The result of this recalculation is the following: 



Fe = 30.86 • km • i? ■ (kpc)"\ /ae + ^T^x^^f" V{>irQ)U{xR) (57) 
and (r R) 

Ve = 30.86 • 10^^ • km • [aqr^ - 3i?S] • (kpc)"^ . (58) 

Both equation are the basis for our numerical investigations, particularly for 
plotting the rotation curves and the spirals, that will be done later. 
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6.2 Spirality of the motion of the test body in the exterior 

The spiral form of the orbiting test body results by superposition of the radial 
motion, described by the radial velocity (jH), and the azimuthal velocity (p7|). Let 
us introduce the notion "spirality" by the definition 

Sp= — =ta,nx, (59) 

Ve 

where x will be called spirality angle. Then the final formula for the spirality reads 



Sp = ta.nx= 7 r ■ (60) 

From this expression we learn that the reversion point of the radial motion is de- 
termined by the quantity 

Rrev = • (61) 

Remembering the relation 

a) Ve^ R and b) V,, = i?$ (62) 

in polar coordinates, we arrive at the following relationship between the spiral angle 
X and the polar angle $ in form of the differential equation 

d$ 

tanx = i?^, (63) 

where tan x is given by (^^ . 

For our further calculations it is convenient to introduce the quantity u = xR. 
Then we arrive at the differential equation 



, d$ 1 / AqY. QiU{u) 

a) — = -^\/as + ^^+ 3 ' with 

b) Qi = 3 , c) P = Aq-3uS, d) Aq ^ aqxr^^ . 

Integration leads to the implicite form $ = $(i?) of the spiral equation ($0 constant 
of integration, mq — xtq) : 



as 



The intervals of the variables are given by 

a) ro ^ i? ^ Rrev , 

(66) 

h) Uq ^ U ^ xRrev ■ 

For a better astrophysical understanding of the spiral function ( |65| ) and particu- 
larly for plotting we are forced to approximate this integral representation by series 
expansion about the expansion center uq according to — uq + rj. Expanding up 
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to the second order in 77, after a rather lengthy calculation we arrive at the result 
« 1) 

Mo 



N 2 

U-UO+ — {u- Uo) 



$ - $0 = ^0 
where Aq and N are defined by 



(67) 



,^ 4 1 . I QiU{uo) 

a) Ao = -f^\as' 



Po V "0 



^0 



3S -AqY.uI + QiuI cos uo~3QiU{uo) (68) 



^ Po ' 2uQ[asul + AQj:ul + QiU{uo)] 

c) Po - Aq-3uo^. 
Solving the quadratic equation (|67|) we arrive at 



p±^p^ - Co + - $0) (69) 

with 

{Nuo - 2)uo 



a) Co 



c) p 



Nuq - 1 



This is the explicite form of the spiral equation received by approximation up to 
the second order in $. 

By further power series expansion we get in first order (linearized form) the 
Archimedes spiral 

a) i? = ro ± S'o(<I> - $0) with 

1 , (71) 

b) ^0 = i^D./p^~Co. 

Hence we find for the distance (at a fixed angle) between two spiral arms the formula 
Ai? = 27r|5ol. (72) 

7 Numerical evaluation of the theory (model Milky 
Way) 

Our next aim is to evaluate numerically for an orbiting test body the general formu- 
las for its radial velocity, azimuthal velocity (rotation curve) and angular velocity as 
well as the formula for its spiral motion (superposition of the radial and azimuthal 
velocities) . 

The formulas presented above show that we need two different categories of 
parameters, namely first the "present" cosmological parameters as adiabatically 
fixed parameters determined by our cosmological model which we treated in our 
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previous papers, and second the individual parameters which partly are rather well 
known by empirical observations. Partly we have to guess two parameters (dark 
matter parameter scalaric bremsheat parameter (3q) in order to be in rough 
accordance with the observations from the motion of the orbiting stars (empirical 
rotation curves). Of course, our simplified model can only reflect some basic features 
of the astrophysics of galaxies. Specifically our further application of the above 
theory refers to a rough model of the Milky Way (Galaxy). 

7.1 Cosmological and individual parameters 

Numerical cosmological values for the "present time" : 

Here we use our previously published results (Schmutzer 2000b): 

CTc = 65.188, 

E = 5.54- 10-19 s, 

(73) 

S = -3.73 • 10-36 ^ 

as = 1.3- 10-35 
Numerical values for rough modelling of the Galaxy: 

Mc = 1.8- 10'''* g (acting central mass , 

rg — 0.8 kpc (radius of the inner central body) , (74) 

>t = 0.04 (kpc)^^ (dark matter parameter) . 

Numerical value of the scalaric bremsheat parameter (analog to sun matter, being 
explained later): 

/3q 4.39 • 10^ . (75) 
Hence the parameter (^l|) and the quantity (|5^) take the form 
a, = 5.89 • 10-30 s-2, 

1 0776 • 10"3 \ ^'^^'^ 
1.302 • 10-5 + "^^^ " kpc j ■ 10-30 s-2 . 

Concluding this subsection let us add some annotations to the scalaric bremsheat 
parameter. According to its definition ( ^6|) it mainly depends on following values 
of the inner central body of the Galaxy: moment of inertia J (p7|), cubic heat 



expansion coefficient ac and heat capacity W {23) which contains the specific heat 
CQ of the central body. The heat consumption factor fs seems to be of the order 
of magnitude 1. Let us further mention that we know the radius tq of the central 
body and roughly its angular velocity lo. Further we realize that the mass of the 
central body cancels. Therefore we are mainly left with the parameters ac and cq . 
Let us for numerical estimates tentatively use the values for the sun, being available 
in astrophysical literature: 



ac™„ « 3.66 • 10-3K \ cg^un « 1-24 • 10« cm^ s-^ K (77) 

Hence we obtain for the sun the value Pq^^^ ~ 4.39 • 103. This is our explanation 
for our previous choice of the parameter Pq. 
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7.2 Plotting 

The above considerations led us to the parameters needed for pfotting, i.e. now we 



are able to plot the radial velocity curve ( |45| ) and (jSg), the rotation curve (17) and 
(19), and the angular velocity curve dH) and (|0l). Fig.l and Fig. 2 show the radial 



course of the radial velocity and the azimuthal velocity V$ (abscissa R in kpc, 
ordinates Vji and V$ in km s^^). Fig. 3 presents the radial course of the angular 
velocity u in s~^. 

Plotting of the spiral of the orbiting body is based on the formula ( |69| ) containing 
the parameters p, Cq and D, which after some substitutions are determined by the 
equations (|7^). A rather lengthy numerical calculation leads to the values 

p = 1.33 kpc, 

Co = 1.49(kpc)^ (78) 

D = -0.259. 

Fig. 4 shows the spiral, where the radial distance is given in kpc. 

I would like to express my thanks to Prof. Dr. A. Gorbatsievich (University of 
Minsk) for scientific discussions and to Mr. E. Schmidt (Jena) for technical advice. 
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Figure 1: Radial course of the radial motion of the star 
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Figure 2: Radial course of the azimuthal velocity of the star (rotation curve) 
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Figure 3: Radial course of the angular velocity 
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